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A consistent combination of covariant density functional theory (CDFT) and Landau- 
Migdal Theory of Finite Fermi Systems (TFFS) is presented. Both methods are in princi- 
ple exact, but Landau-Migdal theory cannot describe ground state properties and density 
functional theory does not take into account the energy dependence of the self-energy and 
therefore fails to yield proper single-particle spectra as well as the coupling to complex con- 
figurations in the width of giant resonances. Starting from an energy functional, phonons 
and their vertices are calculated without any further parameters. They form the basis 
of particle-vibrational coupling leading to an energy dependence of the self-energy and an 
induced energy- dependent interaction in the response equation. A subtraction procedure 
avoids double counting. Applications in doubly magic nuclei and in a chain of superfluid 
nuclei show excellent agreement with experimental data. 



I. INTRODUCTION 

The understanding of the structure of nuclei far from stability with extreme isospin is one the 
most exciting challenges of present nuclear physics. New experimental facilities with radioactive 
nuclear beams make it possible to investigate the nuclear chart to the very limits of nuclear binding. 
A wealth of structure phenomena in exotic nuclei have been reported and the next generation of 
radioactive-beam facilities will present new exciting opportunities to study not only the ground 
states but also excitations and spectra of these strongly interacting many-body systems. This 
situation has stimulated considerable new efforts on the theoretical side to understand the dynamics 
of the nuclear many-body problem by microscopic methods. Exact solutions of the non-relativistic 
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Schrodinger equation based on the bare nucleon-nucleon interaction are used to study very light 
nuclei with ^ < 12 by an "ab initio" approach, modern shell-model calculations based on large scale 
diagonalization techniques and truncation schemes show considerable success in situations where 
configuration mixing calculations are possible, i.e. in light nuclei or in nuclei with single magic or 
doubly magic configurations. For the large majority of nuclei, however, a quantitative microscopic 
description is, so far, only possible by density functional theory (DFT) and its extensions. Although 
DFT can, in principle, provide an exact description of many-body problems{l|, if the exact density 
functional is known, in nuclear physics one is far from a microscopic derivation of this functional. In 
addition, nuclei are self-bound systems. As a consequence of translational invariance the density in 
the laboratory frame is constant in space. Density functional theory in nuclei is therefore based on 
the intrinsic density, a concept that requires additional approximations . The most successful 
schemes of DFT in nuclei use a phenomenological ansatz incorporating as many symmetries of the 
system as possible and adjusting the parameters of the functional to ground state properties of a 
few characteristic nuclei on the nuclear chart. Considerable progress has been reported recently in 
constructing such functionals. For a recent review see J]. 

One of the underlying symmetries of QCD is Lorentz invariance and therefore covariant density 
functionals are of particular interest in nuclear physics. This symmetry not only allows to 

describe the spin-orbit coupling, which has an essential influence on the underlying shell structure, 
in a consistent way, but it also puts stringent restrictions on the number of parameters in the 
corresponding functionals without reducing the quality of the agreement with experimental data. 

Most of the nuclei are superfluid systems and therefore the inclusion of pairing correlations 
is essential for a correct description of structure phenomena in open-shell nuclei [3, B] • Hartree- 
Bogoliubov theory provides a unified description of ph- and pp-correlations on a mean-field level 
by using two type of densities, the normal density matrix p = (a^a) and antisymmetric pairing 

can be combined to the generalized 
density matrix 1Z of double dimension [2]. CDFT theory for superfluid systems is therefore based 
on a generalized Relativistic Hartree-Bogoliubov (RHB) energy density EjiHBiT^]- The same is 
true for the Landau-Migdal theory and for all the methods discussed in this paper. For simplicity, 
however, we restrict all our considerations in this article to the case without pairing correlations, 
^airing correlations can be included on all steps by using super-matrices. Details are given in Ref. 



Only in the applications we present also calculations in isotopic chains of open shell nuclei, 
that include pairing correlations. 

A very successful example of a covariant density functional theory is the Relativistic Hartree- 
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Bogoliubov mode l 11211 . It combines a density dependence through a non-hnear couphng between 
the meson fields [l3] with pairing correlations based on an effective interaction of finite range. 
A large variety of nuclear phenomena have been described over the years within this model: the 
equation of state in symmetric and asymmetric nuclear matter, ground state properties of finite 
spherical and deformed nuclei all over the periodic table 



14( 1 from light nuclei 



1511 to super-heavy 



elements [16], from the neutron drip line, where halo phenomena are observed 



drip line 



13] to the proton 



18| with nuclei unstable against the emission of protons [19]. 



In principle density functional theory can be used for the description of all properties depending 
on the single-particle density. It is therefore not only limited to the description of the ground state 
properties. The same density functionals have also been applied for a very successful description of 
excited stateSj_ such as rotational bands in normal and super-deformed nuclei 
vibrations 



211 ] and collective 



221 ]. Rotations are treated in the cranking approximation providing a quasi-static 
description of the nuclear dynamics in a rotating frame and for the description of vibrations a 
time-dependent mean field approximation is used by assuming independent particle motion in time- 



dependent average fields |23l ]. In the small amplitude limit one obtains the relativistic Random 
Phase Approximation (RRPA) 2^ and in superfluid nuclei the relativistic Quasiparticle Random 
Phase Approximation (RQRPA) [25]. This method provides a natural framework to investigate 
collective and non-collective excitations of ph- (or 2qp) character. It is successful in particular for 
the understanding of the position of giant resonances and spin- or/and isospin-excitations as the 
Gamov Teller Resonance (GTR) or the Isobaric Analog Resonance (lAR). Recently it has been 
also used for a theoretical interpretation of low lying El-strengths observed in neutron rich isotopes 



(pygmy modes) [25l ] and for low-lying collective quadrupole excitations 26|] 



Density functional theory in nuclei is based on intrinsic densities and on the mean field approach. 
Therefore it cannot provide an exact treatment of the full nuclear dynamics. It breaks down in 
transitional nuclei, where the intrinsic density is not well defined and where one has to include 
correlations going beyond the mean field approximation by treating quantum fiuctuations through 
a superposition of several mean field solutions, as for instance in the Generator Coordinate Method 
(GCM) [eI]. But it also provides a poor approximation for the single-particle spectra particularly 
in ideal shell model nuclei such as 208p|-, -^[H^ closed protons and neutron shells. One finds in 
self-consistent mean field calculations usually a considerably enhanced Hartree-Fock gap in the 
single-particle spectrum and a reduced level density at the Fermi surface as compared with the 
experiment. It is well known that this fact is connected with the relatively small effective mass 
in such models. Mahaux and collaborators [2^] have shown that the effective mass in nuclear 
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matter is roughly m* /m ~ 0.8. In finite nuclei it should be modified by the coupling of the single- 
particle motion to low- lying collective surface vibrations. This leads, in the vicinity of the Fermi 
surface, to an enhancement of m* /m ~ 1. Non-self-consistent models with the bare mass ( m* /m 
~ 1) show indeed a single-particle spectrum with a level density close to the experiment. With a 
few exceptions, where the quadrupole motion has been studied within the relativistic Generator 



Coordinate Method (GCM) 23, |29|, applications of covariant density functional theory to the 



description of excited states are limited to relativistic RPA, i.e. to configurations of Ipl/i-nature. 
None of these methods, however, can be applied to the investigation of the coupling to more 
complicated configurations, as it occurs for instance in the damping phenomena causing the width 
of giant resonances. 

Already before density functional theory has been introduced in the sixties for the description 
of quantum mechanical many-body problems by Kohn and Sham [1] Landau has developed in the 
fifties his Fermi Liquid Theory (FLT) [30,] for infinite systems. It has been extended to the Theory 
of Finite Fermi Systems (TFFS) 



3l| by Migdal. This theory provides another very successful 
method for the description of low-lying nuclear excitations [32]. It has several general properties 
in common with density functional theory. First, both theories are know to be exact, at least 
in principle, but in practice, in nuclear physics, the parameters entering these theories have to 
be determined in a phenomenological way by adjustment to experimental data. Second, both 
theories are based on a single-particle concept. DFT uses the mean field concept with Slater 
determinants in an effective single-particle potential as a vehicle to introduce shell-effects in the 

33]. Fermi liquid theory is based on 



exact density functional introduced by Hohenberg and Kohn 
the concept of quasi-particles obeying a Dyson equation, which are defined as the basic excitations 
of the neighboring system with odd particle number. Third, in practical applications both theories 
describe in the simplest versions the nuclear excitations in the RPA approximation, i.e. by a linear 
combination of p/i-configurations in an average nuclear potential. 

However, there are also essential differences between these two concepts. First, in contrast to 
DFT, TFFS does not attempt to calculate the ground state properties of the many-body system, 
but it describes the nuclear excitations in terms of Landau quasi-particles and their interaction. 
Therefore the experimental data used to fix the phenomenological parameters of the theories are 
bulk properties of the ground state in the case of density functional theory, and properties of 
single-particle excitations and of the collective excitations in the case of finite Fermi systems the- 
ory. Second, in DFT the mean field is determined in a self-consistent way and therefore the RPA 
spectrum contains Goldstone modes at zero energy. This is usually not the case in TFFS calcula- 



5 



tions, which are based on a non-relativistic shell-model potential, whose parameters are adjusted 
to the experimental single-particle spectra. Therefore, apart from a few approximate attempts to 
treat the Goldstone modes by adjusting additional parameters in the effective quasiparticle inter- 
action, there is no self-consistency in the RPA calculations of TFFS and the Goldstone modes do 
not separate from the other modes. They are distributed among the low- lying excitations. Third, 
modern versions of TFFS go much beyond the mean field approximation. The coupling between the 
particles and the phonons is investigated with Green's function techniques. Based on the phonons 
calculated in the framework of the RPA one has included particle^phonon coupling vertices and an 



341, 



351] . This leads also to an induced 



energy dependence of the self-energy in the Dyson equation 
interaction in the Bethe-Salpeter equation caused by the exchange of phonons which also depends 
on the energy. The coupling of particles and phonons has also been derived from Nuclear Field 



Theory (NFT) and its extensions [36h38I] . Many aspects of the coupling between the quasi-particles 



and the collective vibrations have been investigated with these techniques 
other kinds of approaches beyond RPA {49, 



39l448l| as well as with 



50| over the years. 
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51 



53l | to find a combination of the ba- 



We give here an overview over recent attempts 
sic ideas of covariant density functional theory and Landau-Migdal theory and show as examples 
corresponding applications. The concept is similar to earlier work in Refs. 
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5a], where specific 



non-relativistic energy functionals have been used to construct a Self-Consistent Theory of Finite 
Fermi Systems. The starting point is a covariant density functional E[p\ widely used in the litera- 
ture. It is adjusted to ground state properties of characteristic nuclei and, without any additional 
parameter, it provides the necessary input of finite Fermi systems theory, such as the mean field 
and the single-particle spectrum as well as an effective interaction between the p/i-configurations 
in terms of the second derivative of the same energy E[p\ with respect to the density. Thus the 
phenomenological input of Landau-Migdal theory is replaced by the results of density functional 
theory. The same interaction is used to calculate the vertices for particle-vibration coupling |5ll |. 



In a second step techniques of Landau-Migdal theory and its modern extensions are used to de- 
scribe the coupling of one- and two-quasiparticle configurations. The main assumption of the 
quasiparticle-phonon coupling model is that two types of elementary excitations - two-quasiparticle 
and vibrational modes - are coupled in such a way that configurations of Iplh^phonon type with 
low- lying phonons strongly compete with simple Iplh configurations close in energy or, in other 
words, that quasiparticles can emit and absorb phonons with rather high probabilities. In this 
way a fully consistent description of the many-body dynamics is obtained. As a result an induced 
additional interaction between single-particle and vibrational excitations provides a strong frag- 
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mentation of the pure RQRPA states causing the spreading width of giant resonances and the 
redistribution of the pygmy strength. 

Two essential approximations are used in this context: First, the Time-Blocking Approximation 
(TEA) 



44 ] ■ that has been extended to systems with pairing correlations (QTBA) in Ref. 
blocks in a special time-projection technique the Ipl/i-propagation through states which have a 
more complex structure than Iplh^phonon. The nuclear response can then be explicitly calculated 
on the lplh+lplh(i^phonon level by summation of an infinite series of Feynman's diagrams. Second, 
a special subtraction technique guarantees, that there is no double counting between the additional 
correlations introduced by particle-vibration coupling and the ground state correlations already 
taken into account in the phenomenological density functional. These two tools are essential for 
the success of this method. TEA introduces a consistent truncation scheme into the Eethe-Salpeter 
equation and without it it would be hard to solve the equations explicitly. The subtraction method 
is the essential tool to connect density functional theory so far used only on the level of mean 
field theory, i.e. on the RPA level, with the extended Landau-Migdal theory, where complex 
configurations are included through particle-vibration coupling. 

The structure of the paper is the following: In Section II we discuss shortly the general formal- 
ism of covariant density functional theory, we introduce in Section III the concept of the energy- 
dependent self-energy S(e) and the vertices of particle-vibration coupling in the relativistic frame- 
work, and we discuss in Section IV the response formalism, the time blocking approximation and 
the subtraction mechanism for the response function. In Section V we present recent numerical 
applications for the calculation of level densities at the Fermi surface and the spreading width of 
the several nuclei. Section VI contains a brief summary and an outlook for future applications. 



II. THE RELATIVISTIC ENERGY DENSITY FUNCTIONAL 



Covariant density functional theory uses the Walecka model [57| as a Lorentz invariant frame- 
work for the formulation of the density functional. In this model the nucleus is described as a 
system of Dirac nucleons coupled to the exchange mesons and the electromagnetic field through an 
effective Lagrangian. The isoscalar scalar cj-meson, the isoscalar vector w-meson, and the isovector 
vector p-meson build the minimal set of meson fields that together with the electromagnetic fie^ 
is necessary for a quantitative description of bulk and single-particle nuclear properties 
The model is defined by the Lagrangian density 

£ = £-N + + C,int- (1) 
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£ N denotes the Lagrangian of the free nucleon 

CN = ^P{i7^d^-m)xP, (2) 

where m is the bare nucleon mass and ip denotes the Dirac spinor. Cm is the Lagrangian of the 
free meson fields and the electromagnetic field 

- \R,a^'' + \mlp,jf - \f,,F^'^, (3) 

with the corresponding masses mo-, mj^, nip, and il^j^, R^u, Ffiu are field tensors (arrows denote 
isovectors and boldface symbols are for vectors in ordinary space). The minimal set of interaction 
terms is contained in Cint 

Cint = -V^r<,aV - - i^r^pppi; - v^r^^^V^- (4) 

with the vertices 

r„ = g^, ^^^ = 5a.7^ f;^ = 5pr7^ r^e=<ir, (5) 

with the coupling constants Qa-, Qui, 9p and (7 (e or for protons or neutrons). Already in the 
earliest applications of the RMF framework it was realized, however, that this simple model with 
interaction terms only linear in the meson fields, does not provide a quantitative description of 
complex nuclear systems. An effective density dependence was introduced [3] by replacing the 
quadratic cj-potential ^m'^a'^ with a non-linear meson coupling potential U{a), which contains 
additional parameters. This particular form of the non-linear potential has become standard in 
applications of RMF models, although additional non-linear interaction terms, both in the isoscalar 
and isovector channels, have also been considered 
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M\ 



^From the model Lagrangian density the classical variation principle leads to the equations of 
motion: 

hf'{id^ + Vp) + m + S]iP = 0. (6) 

If one neglects retardation effects for the meson fields, which is well justified because of the large 
meson masses, a self-consistent solution is obtained when the time-dependent mean-field potentials 



S = gad , = g^uj^ + gpTp^ + qA^ , 
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are calculated at each step in time by the solution of the static Klein-Gordon equations 

- ^(I)m + = ± (V^r^V') , (7) 

where the (+) sign is for vector fields and the (— ) sign for the scalar field. The index m denotes 
mesons and the photon, i.e. (jim = {a , uj^ , f?^ , A^} , and U'{4>m) is derivative of the corresponding 
potential with respect to the meson field. 

In applications to nuclear matter and finite nuclei, the relativistic models are used in the no- 
sea approximation, i.e. the Dirac sea of states with negative energies does not contribute to the 
densities and currents and one uses 

A 

(Vir^V) = (r, t) , (8) 

1=1 

where the sum runs only over the occupied states in the Fermi sea, i.e. vacuum polarization 
effects are neglected. In fact, many effects that go beyond the classical mean- field level are appar- 
ently neglected in the this models: Fock terms, vacuum polarization effects, and the short range 
Brueckner-type correlations. The experimental data to which the meson-nucleon couplings are 
adjusted, however, contain all these effects and much more. It follows that effects beyond the 
mean-field level are implicitly included in the RMF approach by adjusting the model parameters 
to reproduce a selected empirical data set. Vacuum effects, chiral symmetry, nucleon substructure, 
exchange terms, long- and short-range correlation effects are, therefore, effectively included in this 
approach although neither of them can be accessed separately. 

The set of coupled equations ([6]) and d?]) define the relativistic mean field (RMF) model. In 
the stationary case they reduce to a nonlinear eigenvalue problem and in the time-dependent case 
they describe the nonlinear propagation of the Dirac spinors in time [23]. 



RMF models can be also formulated without explicitly including mesonic degrees of freedom. 
Meson-exchange interactions can be replaced by local four-point interactions between nucleons. 
It has been shown that the relativistic point-coupling models |65l467l| are completely equivalent 



to the standard meson-exchange approach. In order to describe properties of finite nuclei on a 
quantitative level, the point-coupling models include also some higher order interaction terms. For 
instance, six-nucleon vertices (tptp)'^, and eight-nucleon vertices (ipip)^ and [('07^'0)(^7'^'(/')]^- 

These relatively simple models turn out to provide a very successful phenomenological descrip- 
tion of the nuclear many-body system all over the periodic table. Relatively few parameters are 
adjusted to ground state properties of a few finite nuclei. At a first glance it is not easy to see 
how can such a simple approach can be so successful. This can be only understood if one considers 
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that this model represents an approximate implementation of Kohn-Sham density functional the- 



ory (DFT) 



111, 3: 



69( 1 . which is successfully employed in the treatment of the quantum many-body 



problem in atomic, molecular and condensed matter physics. 

It is evident that equations of motion Q and ([7]) can also be directly derived form a density 
functional. Using the definition of the relativistic single-nucleon density matrix 

A 

p{v,v\t) = Y,\U^,t)){^,{v',t)\ , (9) 

i=l 

the total energy can be written as a functional of the density matrix p and of the meson fields 

'1 



Ermf [p, ' 



Tr [(ap + /3m) 



d'r + Tr [{Tm<Pm)p] ■ (10) 



The trace operation involves a sum over the Dirac indices and an integral in coordinate space. The 
index m is used as generic notation for all mesons and the photon. 



III. THE ENERGY DEPENDENCE OF THE SELF-ENERGY 

In a relativistic many-body theory the motion of single nucleons in the nuclear medium is 
described by the Dyson equation 

=0, (11) 

where the self-energy is given by 

m* =m + 'Es (12) 

with the scalar part of the self-energy and where the generalized four- vector momentum operator 
has the form 

P^=p^-S^= (i--So,iV + 5]) (13) 
with the vector part E'^ of the self-energy 

S^ = (S°,5]). (14) 

The index 's' in the Eq. (|l'2p denotes that the effective mass is described by the scalar c-meson 
field. In order to characterize ground state properties the stationary Dirac equation has to be 
solved: 

(a(p-S) + /3m* + So)|V) = e|V')- (15) 
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In the general case the full self-energy is non-local in space and also in time. This non-locality means 
that its Fourier transform has both momentum and energy dependence. We therefore decompose 
the total self-energy in a stationary local part and an energy-dependent non-local term: 

S(r, r'; u) = t{v)5{v - r') + E^(r, r'; w), (16) 

where all the components of the self-energy are involved: 

and the index "e" indicates the energy dependence. 

The energy-independent parts of the self-energy correspond to the average fields of the Walecka 
model: 

S,(r) = 5(r), S^(r) = y^(r) (17) 

These fields satisfy the inhomogeneous Klein-Gordon equations, where the sources are determined 
by the respective density and current distributions in a system of A nucleons. 

We assume time-reversal symmetry that means the absence of currents in the nucleus and, thus, 
we find vanishing space-like components of S. The equation of the one-nucleon motion has the 
form: 

{hD + PK{^) + T.l{e))\^)=e\,l:) (18) 

where is the Dirac Hamiltonian 

/id = ap + /3(m + S,) + So (19) 

or, in the language of Green's functions 

{e-hD-liK{^)-T.l{e))G{e) = l. (20) 

It turns out to be useful to work in the shell- model Dirac basis {iV'fc)} which diagonalizes the 
energy-independent part of the Dirac equation: 

/idIV'/c) = efclV'fc)- (21) 
In this basis one can rewrite Eq. (j20p as follows: 

- ek)5ki - ^%{e)]Gik'{e) = Skk', (22) 




FIG. 1: The particle Sp,^,, and the hole components of the relativistic self-energy in the graphical 

representation. Solid and dashed lines with arrows denote one-body propagators for particle (p), hole {h), 
and antiparticle (a) states. Wavy lines denote phonon (/i) propagators, empty circles are the particle-phonon 
coupling vertices 7^ in Eq. The time direction is from the left to the right. 

where the letter indices k,k',l denote full sets of the spherical quantum numbers. 

Obviously, on this stage one needs some model assumptions. The particle-phonon coupling 
model 36( provides a rather simple approximation to describe the energy dependence of S'^(e). 
Within this model 5]'^(e) is a convolution of the particle-phonon coupling amplitude T and the 
exact single-particle Green's function [7Q |: 
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Ki{e) = yl / —Tkvik'{u:)Gk'v{e + u,), (23) 



k'l _ 

where the amplitude 

Tki,i^, (u:) = -y( _JKhjh ^kk'^^i' ) (24) 

is represented in terms of phonon vertexes 7^^ and their frequencies i}^. They are determined by 
the following relation: 

"^ki = X] ^kvik'^Pk'v- (25) 

k'V 

5p is the transition density and in the linearized version of the model V^i'ik' denotes the relativistic 
matrix element of the static residual interaction, i.e. the second derivative of the energy functional 
with respect to the density matrix 

V = ^^^^ (26) 
opop 
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In the linear approximation 5p is not influenced by the particle-phonon coupling and can be 
computed within relativistic RPA or QRPA. The linearized version implies also that the energy- 
dependent part of the self-energy (123p contains the mean field Green's function G{s) = [e — 
instead of the exact Green's function G. The graphical representation of the self-energy is given 
in Fig. [TJ In contrast to the non-relativistic case, where one has occupied states below the Fermi 
surface (hole states h) and empty states above the Fermi surface (particle states p) we now have 
according to the no-sea approximation in addition empty states with negative energies in the Dirac 
sea (anti-particle states a). Particle and hole components are drawn assuming all the possible 
types of intermediate states. 

Eq. (|22p contains off-diagonal elements of the matrix with relatively large ener gy d enomi- 



3|] that it 



nators. It has been shown by explicit calculations within the non-relativistic approach 
is justified to use the diagonal approximation: 

Ellis) = 5uiK{e). (27) 



Thus, within the diagonal approximation of the self-energy (j27p the exact Green's function G is 
also diagonal in the Dirac basis and the Dyson equation forms for each k a non-linear eigenvalue 
equation 

{e-ek-n{e))Gk{e) = l. (28) 
The poles of the Green's function Gk{e) correspond to the zeros of the function 

f{e)=e-ek-n{e). (29) 



In Refs. 



511 ] it is shown how this problem can be solved by a matrix diagonalization. For each 



quantum number k there exist several solutions e'^'^ characterized by the index A. Because of the 
coupling to the collective vibrations the single-particle state k is fragmented as it will be shown in 
the application in section IV A[ 



IV. THE RESPONSE FUNCTION 

Nuclear dynamics of an even-even nucleus in a weak external field is described by the linear 
response function i?(14, 23), where 1 = {ki,ti} combines the quantum numbers k and the time. 
This response function is the solution of the Bethe-Salpeter equation (BSE) in the ph channel: 

i?(14, 23) = i?°(14, 23) - i ^ i?°(16, 25)^(58, 67)i?(74, 83), (30) 

5678 
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with the free response i?''(14, 23) = 3)(j'(4, 2), where function G is the exact single-particle 
Green's function. The summation over the number indices 1, 2, . . . implies an integration over the 
respective time variables, and V is the effective interaction irreducible in the p/i-channel. This 
interaction is determined as a variational derivative of the full self-energy S with respect to the 
exact single-particle Green's function: 

V-(14,23,..||M. ,31, 



Since the self-energy in Eq. (116p has two parts S = S -|- S^, the effective interaction V in Eq. (j30l 
is a sum of the static RMF interaction V (j26p and time-dependent terms 



After a Fourier transformation in time, this time dependence leads to an energy-dependent inter- 
action V^. In the Dirac basis (1211) it has the form: 



V^,i^,{u:,e,e') = V ^2V^-Jl_^. (33) 

lJ,,a ^ " 

where a = +1 for empty states and —1 for occupied states (for details see Ref. [s^). The Bethe- 
Salpeter equation (j30p contains the exact Greens' function G. In order to simplify this equation 
for the further analysis the G is expressed it in terms of the mean field Green's function G. From 
Eq. (j22p we derive the Nambu form for it: 

(5-1(1,2) = G-i(l,2) + S^(l,2), (34) 

which reads in Fourier space as 

G'fcifc2(e) = . (35) 

Introducing ^°(14,23) = (5(1, 3)G'(4, 2) one can rewrite Eq. ([30]) as follows: 

R = R° - iR^WR, (36) 
where is a new interaction of the form 

W = V + W^ (37) 

with 

^^^^(14, 23) = ^^^(14,23) +iS^(l,3)(5-H4,2) +i(5-i(l,3)S"(4,2) -iS'^(l,3)S^(4,2). (38) 
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FIG. 2: Bethe-Salpeter equation for the p/i-response function R in graphical representation. Details are 
given in Fig. [T]and the small black circle means the static part of the residual p/i-inter action 



The graphical representation of the Eq. ()36p is shown in Fig. [2j 

In addition to the static interaction V the effective interaction W contains diagrams with energy- 
dependent self-energies and an energy-dependent induced interaction, where a phonon is exchanged 
between the particle and the hole. As discussed in Refs. 



44 



46 



52] the term zS''(3, 1)S^(2,4) has 



to neglected in the time blocking approximation if one neglects backward-going propagators caused 
by the particle-phonon coupling. This is a reasonable approximation applied and discussed in many 



non-relativistic models (see e.g. Refs. [39l448l| and references therein). We have to emphasize, 
however, that all the RPA ground state correlations are taken into account, because it is well 
known that they play a central role for the conservation of currents and sum rules. 

It turns out that both the solution Eq. (136p R and its kernel W are singular. Another difficulty 
arises because Eq. (|36p contains integrations over all time points in the intermediate states. This 
means that many configurations are contained in the exact response function which are actually 
more complex than Iplh phonon. Tselyaev has introduced in the Ref. |4j] the Time Blocking 
Approximation (TBA), a special time-projection technique to block the p/i-propagation through 
these complex intermediate states. In this way one obtains after a Fourier transformation in time 
a relatively simple algebraic equation: 

R{uj) = R^{uj) + R^{uj)W{uj)R{uj), (39) 

where 

Wk^kiMkai^) = ykikiMk-i + ^kikiMk-ii^) - ^kikiMkai^) (40) 

and 



RkikiMkii^) = Rkik2{^)Skik3Sk2k4- 



(41) 
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Rkik2{^) is the mean field propagator: 



Rph{oj) = ^ , Rah{uj) = ^ , (42) 

£ph — W Eah — W 
Rhp{^) = Rhaiuj) = (43) 



^ph = ^p ^ and ^ is the particle-phonon couphng amphtude with the following components: 



'^h"hyh"h' 



fi /I* /i fi* 

f 1pp"^p>p" sr^ 

+ Ohh' I 2^ ^ pr- + 2^ 



p" ^ a" 

"Ipp'^hh' _|_ ^p'p^h'h 



-UJ — Ep' + Eh — OJ — Ep + Eh' — 



(44) 



M A** 
^h"h'yh"h' 



iO -Ea+ Eh" - 

W - + efc - ^LO- Ep" +Eh-ni' 

a" p" ' 

fj, fi* 11*11 
^aa'^hh' _|_ ^a'a^h'h 



.OJ - Ea' + Sh-^f' ' LO -Ea+Eh' -^f'^y 

As in Fig. [T]the indices p,a and h denote the particles, antiparticles and holes in the Dirac basis. 



The amplitudes ^ph',ha, ^ah',hp are neglected, because they have only are small effect (see Ref. |5l|). 



The amplitudes ^pp'^hh' and ^hh',pp' are also disregarded within this approximation. Therefore, 
ground state correlations are taken into account only on the RPA level due to the presence of the 
ypp',hh', yhh',pp' terms of the static interaction in the Eq. (p9]) . By definition, the response function 
R[ijj) in Eq. (j39p contains only configurations which are not more complex than Iplh phonon. 

In Eq. ()39p $(0) is subtracted from ^(oj). This corresponds to the subtraction procedure devel- 
oped by Tselyaev in Ref. [47]. It considers the fact that the effective interaction V being adjusted 
to experimental data of the ground state contains effectively many correlations and, in particular, 
also admixtures of phonons at the energy u = 0. In the present method, all correlations entering 
through the admixture of phonons are taken care of by the additional interaction term ^{lo). To 
avoid double counting in the effective interaction the part <J>(0) is therefore subtracted. This means 
only the energy dependence of the phonon coupling is effectively taken into account. 

To describe the observed spectrum of the excited nucleus in a weak external field D, as for 
instance a dipole field, one needs to calculate the strength function: 

S{E) = -- lim Im UDD{E + iA), (46) 
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expressed through the polarizabihty 

UDDii^) = D^R{io)D (47) 

A finite imaginary part A of the energy variable is introduced in the calculations for convenience 
in order to obtain a more smoothed envelope of the spectrum. This parameter has the meaning 
of an additional artificial width for each excitation. This width emulates effectively contributions 
from configurations which are not taken into account explicitly in this approach. 

In order to calculate the strength function it is convenient to convolute Eq. (j39p with an external 
field operator and introduce the transition density matrix 6p in the external field D: 

{^)Dk,k,, (48) 

Using Eq. (p9]l we find that 5p{uj) obeys the equation 

6p{uj) = 6p^{io) + ^°(a;) (v + $(w) - $(0)) <5p(w), (50) 
and the strength function is expressed as 

S{E) = ^im^ImTr[D'<Sp{E + iA)]. (51) 

V. APPLICATIONS 



711 ] is used for the 



For the following applications discussed in this section the parameter set NL3 
covariant energy functional. For superfiuid nuclei we use in the pairing channel a simple monopole 
force with the strength parameters Gr (r = p, n) adjusted to experimental gap parameters for 
protons and neutrons. The cut-off energy in the pairing channel is 20 MeV both for protons as 
well as for neutrons. The parameter set NL3 has been adjusted to ground state properties of a 
few spherical nuclei more than ten years ago. In numerous applications it has been shown that 
it provides on the mean field level a very good description of ground states and excited states 
of nuclei all over the periodic table In a recent investigation {t^ its parameters have been 

slightly modified and several small deficiencies have been eliminated. 

A. Single-particle spectra in the Pb-region 

In this section we discuss the changes of the single-particle spectra of the odd mass nuclei 207pb, 
^°^Pb, 2'^'^Tl and ^''^Bi if the coupling to low lying collective vibrations of the surface is taken into 
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TABLE I: Energies e). and spectroscopic factors of the dominant neutron levels in ^"^Pb calculated 
in the strongly restricted particle-phonon space, pha denotes full the calculation, pa {h) is the version 
without backwards going terms, and ph is the version without contribution of the antiparticle states in the 
self-energy (see text for details). 



State k 


efc, MeV 




eP, MeV 










Particle 




pha 


pa 


ph 


pha 


pa 


ph 


2g9/2 


-2.50 


-2.85 


-3.14 


-2.88 


0.89 


0.92 


0.89 


lill/2 


-2.97 


-2.82 


-3.20 


-2.90 


0.94 


0.97 


0.94 


ljl5/2 


-0.48 


-1.16 


-1.33 


-1.21 


0.70 


0.74 


0.70 


3d5/2 


-0.63 


-0.96 


-1.05 


-0.98 


0.93 


0.94 


0.93 


4sl/2 


-0.36 


-0.88 


-0.92 


-0.89 


0.93 


0.93 


0.93 


2g7/2 


-0.56 


-0.71 


-0.90 


-0.76 


0.92 


0.94 


0.92 


3d3/2 


-0.02 


-0.35 


-0.42 


-0.37 


0.93 


0.93 


0.93 


Hole 




pha 


h 


ph 


pha 


h 


ph 


3pl/2 


-7.66 


-7.67 


-7.40 


-7.70 


0.96 


0.98 


0.96 


2f5/2 


-9.09 


-8.97 


-8.71 


-9.02 


0.93 


0.96 


0.93 


3p3/2 


-8.40 


-8.20 


-7.87 


-8.22 


0.90 


0.94 


0.90 


lil3/2 


-9.59 


-9.30 


-9.07 


-9.36 


0.90 


0.92 


0.89 


2f7/2 


-11.11 


-10.20 


-9.98 


-10.22 


0.72 


0.76 


0.72 


(lh9/2)i 


-13.38 


-13.32 


-13.23 


-13.34 


0.52 


0.47 


0.53 


(lh9/2)2 




-12.48 


-12.42 


-12.49 


0.31 


0.39 


0.29 



account. In order to keep the numerical effort in reasonable limits in a first investigation only the 
most collective phonons with spin and parity J'^ = 2^", 3^, 4+, 5^, 6^ below the neutron separation 
energy and a reduced number of single-particle states with positive energy (particles or holes) is 
taken into account in the solution of the Dyson equation (j28p . This reduces strongly the number 
of poles in the self-energy of Eq. (j23p . The numerical results obtained in these investigations are 
compiled in the Table [H For the first shell of neutron levels above ('particle') and below ('hole') 
the Fermi level three versions are given: in the version pha the index n in Eq. (|23p includes all 
contributions from intermediate states above the Fermi level p, below the Fermi level h and in the 
Dirac sea a. Version pa (for particles) or h (for holes) excludes the backward going diagrams, and 
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the third version ph does not contain antiparticle intermediate states in (j23p . In this way, one 
can see that the effects of ground state correlations (GSC) caused by the particle-phonon couphng 
and neglected in the second version are significant and it is essential to take them into account 
in a realistic calculation. On the other hand, the contribution of the antiparticle subspace to the 
self-energy is quantitatively not of great importance. This can be understood by the large values 
of the energy denominators in Eq. (I23p for these configurations. Thus it is justified to disregard 
them in the following calculations. Notice, however, that version ph does not eliminate the effects 
of the Dirac sea completely since the phonon vertices still contain this contribution. As it has been 
discussed in Ref. [2j] these terms play an important role in a proper treatment of relativistic RPA. 
Otherwise it is not possible to obtain reasonable properties for the isoscalar modes within RRPA. 

Next we show results where the contribution of the antiparticle subspace to the self-energy are 
neglected. In this case one is able to enlarge the particle-hole basis considerably by taking into 
account particle-hole configurations far away from the Fermi surface. This increases the collectivity 
of the phonons and, consequently, the strength of the particle-vibrational coupling. The phonon 
basis was also enriched by including higher-lying modes up to 35 MeV. Solving the Dyson equation 
one finds a fragmentation of the single-particle states and a corresponding reduction of the single- 
particle strength. For the levels one major shell below and one shell above the Fermi surface one 
finds always one dominant level, which is shifted against the corresponding single-particle energy 
without particle-phonon coupling. Almost all the levels are moving downwards providing thus 
a considerably better agreement with experimental energies then the pure RMF states. In the 
next shells further away from the Fermi surface almost all the single-particle levels turn out to be 
strongly fragmented due to phonon coupling and it is no longer possible to determine the dominant 
levels in these shells, in other words, the concept of Landau quasi-particles is defined only in the 
neighborhood of the Fermi level and it breaks down at larger distances. 

In Fig. [3] we show as an example the 3d3/2 and the lji5/2 levels in the nucleus ^''^Pb. In both 
cases the single-particle strength is distributed over about two thousand states but most of them 
are vanishingly small. Thus only the states with the strength exceeding 10^^ are drown. The state 
3d3/2 has a pronounced single-particle structure with a single-particle strength close to 0.9. On 
the other side the lji5/2 is more fragmented. The experimental strength of the dominant levels are 
shown with dashed lines. 

To illustrate the shifts in the level schemes of the dominant poles as compared to the RMF 
results we show as an example in Figs. H] the single-particle spectrum for neutrons. The spectrum 
calculated with the energy-dependent correction (RMF+PVC) demonstrates a pronounced increase 
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FIG. 3: Single-particle strength distribution for the 8^3/2 (left panel) and lji5/2 (right panel) states in ^°^Pb 
obtained in the calculations (solid lines) and the experimental strengths of the respective dominant levels 
(dashed lines). 
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Neutron states in Pb 



FIG. 4: Neutron single-particle states in Pb^°^: the pure RMF spectrum (left column), the levels computed 
within RMF with allowance for the particle-vibration coupling (center) and the experimental spectrum 
(right). 



of the level density around the Fermi surface of ^"*^Pb compared the pure RMF spectra. In some 
cases the order of levels is inverted and the observed sequence is reproduced as for instance for 

the lji5/2 ^-i^d the 3d^/2 neutron states. Another and more important example is the inversion of 
the 2^(9/2 and liii/2 neutron states which reproduces the spin of the ^'^^Pb ground state. 
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In order to quantify these results we calculate the average distance between two levels in the 
spectrum shown in Fig. H One finds for the neutrons 1.0 (RMF), 0.83 (RMF+PCV) and 0.76 
(EXP) in units of MeV. This corresponds to a level density of 1.0 (RMF), 1.20 (RMF+PCV) and 
1.31 (EXP) in units of MeV~^. The level density in the neighborhood of the Fermi surface is 
therefore in RMF-calculations by a factor 0.76 smaller than the experimental value. Taking into 
account particle- vibrational coupling we find only a reduction of 0.92. Assuming an effective mass 
close to 1 for the experiment, and taking into account that the level density at the Fermi surface 
is proportional to m*/m, this corresponds to an effective mass m* /m ~ 0.76 for the RMF and 
m* /m w 0.92 for the RMF+PCV calculations. For the protons the situation is similar. 



Jaminon and Mahaux have discussed in Refs. 

mm 

the concept of the effective mass in the 
case of RMF theory. On one side one has the well known Dirac mass 



which is determined by the scalar field S^. Since we do not use an isovector scalar field for 
the present parameter set NL3 the Dirac mass is in these calculations identical for protons and 
neutrons. However, this quantity should not be compared with the effective mass determined 
empirically from a non-relativistic analysis of scattering data and of bound states. From a non- 
relativistic approximation of the Dirac equations one finds that the mass 



should be used for this purpose. Here Sq is the time-like component of the Lorentz vector field 
determined by the exchange of uj- and /?- mesons. 

In symmetric nuclear matter we find for NL3: rriDlra = 0.60 and meff/m = 0.67. The latter 
value is smaller then the values m* /m ~ 0.71 for protons and m* /m ~ 0.76 for neutrons deduced 
from the calculated spectrum around the Fermi surface in simple RMF theory. Following similar 
arguments we would obtain for RMF+PVC calculations an average effective mass of 0.89. This is 
obviously still too low as compared to the experimental value. 

On the other hand, around the Fermi surface where relativistic kinematic effects are not sig- 
nificant the RMF+PVC spectrum can be characterized by the effective mass deduced from the 
Schrodinger equation which is a non-relativistic limit of the Dirac equation (jlSp . In this approxi- 
mation one can calculate the state-dependent E-mass fh/m^^'^^ which is the inverted spectroscopic 
factor of the dominant level A: 




niD = m + Ss(r) 



(52) 



rueff = m-T,o 



(53) 
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For the calculated RMF+PVC spectrum the averaged E-masses are 1.26 for neutrons and 1.41 for 
protons if one takes into account all the states with spectroscopic factors larger then 0.5, i. e. good 
single-particle states. Thus, the energy dependence of the self-energy increases the RMF neutron 
and proton effective masses up to the values 0.96 and 1.0, respectively. 

Although the problem of particle-vibration coupling in nuclei has a long history and it was 
considered in a number of works, most of them are based on a non-relativistic treatment of the 
nuclear many-body problem. Only in a relatively recent investigation in Ref. [3] a correction of 
the RMF single-particle spectrum was undertaken in a phenomenological way assuming a linear 
dependence of the self-energy near the Fermi surface. The corresponding coupling constants were 
determined by a fit to nuclear ground state properties. Despite the fact that the present approach is 
fully microscopic without any additional parameter adjusted to experiment it shows good agreement 
with the results of Ref. 7a] for the spectrum of ^'^'^Pb. The shift caused by the phenomenological 



particle-vibrational coupling in Ref. 



Non-relativistic microscopic investi] 
two major groups. The first group [34 



is only slightly larger than in the present investigation. 
:ations of particle-vibrational coupling can be divided into 
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7C 



77[ uses a phenomenological single-particle input 



to reproduce the experimental spectrum and has therefore to exclude the contribution of the 
particle-vibration coupling from the full self-energy to find the 'bare' spectrum. Usually these 
older approaches take into consideration only a relatively small number of collective low-lying 



phonons and use a particle-vibration coupling model 



3a| . This restriction to only low-lying modes 



produces shifts less then 1 MeV. However, as it was shown in Ref. 



35[, enlarging of the phonon 



space with high-lying vibrations leads to very strong shifts of the single-particle levels up to 4 MeV, 
and no saturation is observed with respect to the dimension of the phonon space. 

The second group of approaches (see, for instance Refs. QjS]) starts from a self-consistent 
Hartree-Fock description and applies perturbation theory to calculate the particle-vibration con- 
tribution to the full self-energy. In such self-consistent methods it is more justified to enlarge the 



phonon space. It was shown, for instance, in 



□ 



that the contribution of the isovector modes is 
noticeably smaller than the isoscalar ones. The detailed investigation of the relative importance 
of the high multipole states was performed in [78]. Because of the larger phonon space the typical 
shifts of the single-particle levels in 208p]-, about 1-2 MeV. 

As for the spectroscopic factors, all the approaches predict similar values because these factors 
are not very sensitive to the details of the calculation schemes. 
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B. The strength functions of collective excitations in closed shell nuclei 

The solution of the Bethe-Salpeter equation ([39|) ahows to calculate the nuclear response to 
external multipole fields and the strength functions of the corresponding collective excitations. As 
in the last section we show applications based on the density functional NL3 with a monopole force 
in the pairing channel. A small artificial width of 200 keV is introduced as an imaginary part of the 
energy variable to to have a smooth envelope of the calculated curves. The energies and amplitudes 
of the most collective phonon modes with spin and parity 2+, 3^, 4^, 5~, 6"*" are calculated with 
the same restrictions and selected using the same criterion as in the last section and in many other 
non-relativistic investigations in this context. Only the phonons with energies below the neutron 
separation energy enter the phonon space since the contributions of the higher-lying modes are 
found to be small. 

On all three stages of these calculations the same energy functional, i.e. the same relativistic 
nucleon-nucleon interaction V (I26p has been employed. The vertices T^^/.^ ^251) entering the term 
<^(cj) in Eq. ()44p are calculated with the same force. Therefore no further parameters are needed. 
The scheme is fully consistent. 

The subtraction procedure developed by Tselyaev in the Ref. 47]] for the self-consistent scheme 



removes the static contribution of the particle-phonon coupling from the p/i-interaction. It takes 
into account only the additional energy dependence introduced by the dynamics of the system. It 
has been found in the calculations of Refs. 



48| | as well as in the calculations of the Ref. [Ul] that 
within the relatively large energy interval (0 - 30 MeV) the subtraction procedure provides a rather 
small increase of the mean energy of the giant dipole resonance (0.8 MeV for lead region) and gives 
rise to the change by a few percents in the sum rule. This procedure restores the response at zero 
energy and therefore it does not disturb the symmetry properties of the RRPA calculations. The 
zero energy modes connected with the spontaneous symmetry breaking in the mean field solutions, 
as for instance the translational mode in the dipole case, remain at exactly the same position 
after the inclusion of the particle-vibration coupling. In practice, however, because of the limited 
number of oscillator shells in the calculations this state is found already in RRPA without particle- 
vibration coupling at a few hundreds keV above zero. In cases, where the results depend strongly 
on a proper separation of this spurious state, as for instance for investigations of the pygmy dipole 
resonance in neutron rich systems 



53|, 



80(] one has to include a large number of p/i-configurations 
in the RRPA solution. 

In Fig. [5] we show the calculated strength functions for the isoscalar monopole resonance in 
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FIG. 5: Isoscalar monopole resonance in ^'^^Pb and ^^^Sn obtained within two approaches: RRPA (dashed 
Une) and RRPA with particle-phonon couphng RRPA-PC (sohd hne). Both computations have been per- 
formed with relativistic Hartree (RH) mean field and employ NL3 parameter set for RMF forces. 

TABLE II: Lorentz fit parameters of isoscalar EO strength function in ^°^Pb and ^'^^Sn calculated within 
RRPA and RRPA extended by the particle-phonon coupling model (RRPA-PC) as compared to experimental 
data. The fit has been carried out in the interval from i?„ to roughly 20 MeV 







<E> (MeV) 


r (MeV) 




RRPA 


14.16 


1.71 


208pb 


RRPA-PC 


14.05 


2.36 




Exp. [81] 


13.73(20) 


2.58(20) 




RRPA 


16.10 


2.63 


"2Sn 


RRPA-PC 


16.01 


3.09 



^°Pb and ^ Sn. The fragmentation of the resonance caused by the particle-phonon coupling is 
clearly demonstrated although the spreading width of the monopole resonance is not large because 
of a strong cancellation between the self-energy diagrams and diagrams with the phonon exchange 
(see Fig. [2]). This fact has also been discussed in detail in Refs. js^, and it is not disturbed by 
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FIG. 6: The isovector El resonance in "°Pb. Details are given in the text. 

the subtraction procedure because this cancehation takes place as well in ^(u) as in <I>(0). 

In order to compare the spreading of the theoretical strength distributions with experimental 
data we show in Table |TI] mean energies (E) and widths parameters F obtained by fitting the 
theoretical strength distribution in a certain energy interval to a Lorentz curve in the same way 
as it has been done in the experimental investigations. The experimental values shown in Table Ull 
are derived in the Ref. 811] from the evaluation of a series of data obtained in different experiments 
for the isoscalar monopole resonance in ■^'^^Pb. 



Figs. [6] and [7] present calculated photoabsorption cross sections 

9hc 



(55) 



for the isovector dipole resonance in ^''^Pb and in ^'^^Sn. The left panels give the results ob- 
tained within the non-relativistic approach with a Woods-Saxon (WS) single-particle potential 



48l |. They are compared in the right panel with 



and Landau-Migdal (LM) forces described in Ref. 

the relativistic fully consistent theory of Ref. 5^. (R)RPA calculations are shown by the dashed 
curves, (R)RPA extended by the phonon coupling ((R)RPA-PC) calculations - by the thick solid 



curves. In Fig. [6] we have also displayed experimental data with error bars taken from Ref. 82l |. 
In both calculations, relativistic (right panel) and non-relativistic (left panel), the continuum is 
taken into account only in a discrete approximation, which is very reliable for heavy nuclei. As 
discussed in Ref. [14], the Dirac equation ()2ip is solved in an oscillator basis. To make the com- 
parison reasonable calculations within the non-relativistic framework have been performed with 
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FIG. 7: The same as in Fig.[6]but for "^Sn. 

box boundary conditions for the Schrodinger equation in r-space which ensures completeness of 
the single-particle basis. 

The corresponding Lorentz fit parameters in the two energy intervals: Bn — 25 MeV and — 30 
MeV (Bn is the neutron separation energy) are included in Table IIIII and they are compared with 



the data of Ref. 



83| . We notice that the inclusion of particle-phonon coupling in the RRPA 



calculation induces a pronounced fragmentation of the photoabsorption cross sections, and brings 
the width of the GDR in much better agreement with the data, both for ^"^^Pb and ^^^Sn. 

The fragmentation of the resonance introduced by the particle-phonon coupling is clearly demon- 
strated in both cases. Also, one finds more or less the same level of agreement between theory and 
experimental data for these two calculations. In the case of the isovector El resonance in ^^^Sn 
this is, however, not so clear because the cross section and the integral characteristics of the reso- 
nance obtained in the experiment of Ref. 8^ are given with relatively large error bars. In ^"^^Pb 
the self-consistent relativistic approach reproduces the shape of the giant dipole resonance much 
better than the non-relativistic one although the whole resonance is about 0.5 MeV shifted to lower 
energies with respect to the experiment. As one can see from the Fig. [6] and Table HIH we observe 
some shift already in the RRPA calculation, which is determined by the properties of the NL3 



forces. Improvement of the forces, for instance, the use of the density dependent versions 
of the RMF should bring the El mean energy in better agreement with the data. 



84], 



85| 
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TABLE III; Lorentz fit parameters in the two energy intervals: Bn ~ 25 MeV and — 30 MeV, for tlie 
El plioto absorption cross sections in ^°*Pb and ^'^^Sn, calculated with the RRPA, and with the RRPA 
extended to include the particle-phonon coupling (RRPA-PC), compared to data. 

Bn - 25 MeV - 30 MeV 







<E> 


r 


EWSR 


<E> 


F 


EWSR 






(MeV) 


(MeV) 


(%) 


(MeV) 


(MeV) 


(%) 




RRPA 


13.1 


2.4 


121 


12.9 


2.0 


128 


208p|2 


RRPA-PC 


12.9 


4.3 


119 


13.2 


3.0 


128 




Exp. [82] 


13.4 


4.1 


117 






125(8) 




RRPA 


14.7 


3.3 


116 


14.5 


2.6 


126 


i32Sn 


RRPA-PC 


14.4 


4.0 


112 


14.6 


3.2 


126 




Exp. [83] 


16.1(7) 


4.7(2.1) 


125(32) 









However, there is an essential difference between the fully self-consistent relativistic calculations 



48l | one introduces 



and the non-relativistic approach: in non-relativistic approach discussed in Ref. 
on all three stages of the calculation phenomenological parameters, which have to be adjusted to 
experimental data: first, the Woods-Saxon parameters as, for instance, the well depth are varied 
to obtain single-particle levels close to the experimental values, second, one of the parameters of 
the Landau-Migdal force is adjusted to get phonon energies at the experimental positions (for each 
mode) and, third, another Landau-Migdal force parameter is varied to reproduce the centroid of 
the giant resonance. Although the varying of the parameters is performed in relatively narrow 
limits, it is necessary to obtain realistic results. In contrast, in the relativistic fully consistent 
approach no adjustment of additional parameters is necessary. Of course, the underlying energy 
functional has been determined in a phenomenological way by a fit to experimental ground state 
properties of characteristic nuclei. However, it is of universal nature and the same parameters are 
used for investigations of many nuclear properties all over the periodic table. The predictive power 
of this scheme is therefore much higher than that of the present semi-phenomenological approach 



discussed, for instance, in Ref. 



48|. 
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C. Collective excitations in systems with pairing 

Pairing correlations play an essential role in all open shell nuclei and apart from the vicinity 
of the very few doubly magic configurations nuclei show superfluidity all over the periodic table. 
In a theoretical description this fact can be taken into account by Bogoliubov's quasiparticles. 
Many-body theories for normal systems are thus relatively easily extended to the case of superfluid 
nuclei. Combining creation and annihilation operators and a to a two-component operator, 
operators of the type a~^a, a'^a'^ , aa are replaced by super-matrices of rank 2 and the form of the 
equations stays nearly unchanged. This very elegant method has been introduced already half a 



century ago in Ref. 



8a] and over the years it has been used for various many-body approximation 



schemes in non-relativistic systems as for instance in Refs. 



Id, 



87 



89l |. Response theory with the 



Time Blocking Approximation (TEA) introduced for normal non-relativistic systems in Ref. 



has been extended in Refs. 



471, 



Ml 



43] to QTBA for superfluid systems and in Ref. [li] to RQTBA 
for relativistic superfluid systems. 

In the following we discuss several applications of RQRPA and of RQTBA in the chain of spheri- 
cal even-even semi-magic nuclei with Z = 50. We show calculations of the isovector dipole spectrum 
in the giant dipole resonance region and in the low-lying energy region in the two approximations. 
As discussed before the effective interaction <I>(cj)-<I>(0) takes into account only the additional 



energy dependence introduced by the dynamics of the system. It has been found in relativistic 
as well as in non-relativistic calculations [48] that within a relatively large energy interval (0 - 30 
MeV) the subtraction procedure provides a rather small but noticeable increase of the mean energy 
of the giant dipole resonance (about 0.7 MeV for tin region) and gives rise to changes by a few 
percents in the sum rule. The absolute value of the energy shift produced by the subtraction of ^(0) 
in Eq. ()40p is comparable with but not exactly equal to the absolute value of the shift produced 
by the dynamical part of the interaction amplitude ^(co) which always reduces the mean energy 
of the resonance. The subtraction procedure restores the response at zero energy and, therefore, 
it does not disturb the symmetry properties of the RQRPA calculations. The zero energy modes 
connected with the spontaneous symmetry breaking in the mean field solutions, as, for instance, 
the translational mode in the dipole case, remain at exactly the same positions after the inclusion 
of the quasiparticle-vibration coupling. 

In Fig. [8] we show dipole spectra for the tin isotopes ^^^Sn, ^^"Sn, ^'^'^Sn. The right panels 
show the photo absorption cross section (j55[) which is determined by the dipole strength function 



Sei- It is calculated, analogously to Eq. (pT]) . with the usual isovector dipole operator. The left 
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FIG. 8: The calculated dipole spectra for the heavier tin isotopes "'^"'^^Sn, -'^^^Sn, ^^°Sn, compared to data of 
Ref. [9^ for ^^^'^^°Sn. Right panels (b, d, f): photo absorption cross sections computed with the artificial 
width 200 keV. Left panels (a, c, e): the low- lying portions of the corresponding spectra in terms of the 
strength function, calculated with 20 keV smearing. Calculations within the RQRPA are shown by the 
dashed curves, and the RQTBA - by the solid curves. 



panels show the low-lying parts of the corresponding spectrum in terms of the strength function. 
A small imaginary part of 20 keV is used for the energy variable, in order to see the fine structure 
of the spectrum and sometimes individual levels in this region. RQRPA calculations are shown by 
dashed curves and the RQTBA by the solid curves. Experimental data are taken from the EXFOR 
database [9C | . 



These figures clearly demonstrate how the two-quasiparticle states, which are responsible for 
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the spectrum of the RQRPA excitations, are fragmented through the couphng to the collective 
vibrational states. The effect of the particle- vibration coupling on the low- lying dipole strength 
below and around the neutron threshold within the presented approach is shown in the left panels 
of the Fig. [8l Such calculations give us an example how the low-lying strength develops with the 
increase of the neutron excess. It is also found that the presence of pairing correlations causes a 
noticeably stronger fragmentation of both the GDR and the PDR modes as compared to the case 
of a normal system discussed above. This effect has the two reasons. First, pairing correlations 
lead to a diffuseness of the Fermi surface and, thus, increase the number of possible 2qp ^phonon 
configurations, and second, pairing correlations cause a considerable lowering of the energies and 
increased transition probabilities of the lowest 2+ states. In spherical open-shell medium mass 
nuclei the highly collective first states appear at energies around 1 MeV (and they are usually 
well reproduced in RQRPA joi]) whereas in magic nuclei and often in nuclei near the shell closures 
they appear much higher, at about 3-4 MeV and have considerably reduced transition probabilities. 
This causes a strong configuration mixing in the case of presence of very low- lying vibrational states. 
These modes admix to others, in particular, to the GDR and the PDR and the lower their energies 
and the higher their transition probabilities are, the stronger fragmentation they cause. 

A systematic analysis of the transition densities of the RQRPA and the RQTBA states shows 
that the 2qp transition densities in the broad low-lying energy region dominated by the frag- 
mentation of the RQRPA pygmy mode have a very similar behavior as the initial RQRPA state: 
proton and neutron components oscillate in phase in the nuclear interior and neutron components 
dominate on the surface in nuclei with noticeable neutron excess. 



VI. CONCLUSIONS 



We have given an overview over recent efforts to combine two theoretical methods for the 
description of the quantum-mechanical many-body problem of nuclear physics, Covariant Density 
Functional Theory (CDFT) and Landau-Migdal Theory of Finite Fermi Systems (TFFS). Both 
methods are very successful and they are claimed to provide in principle an exact description. In 
practice, however, there are limitations. Both methods use phenomenological input. DFT can only 
be applied to physical quantities, which can be expressed in terms of the single-particle density and 
in self-bound systems such as nuclei DFT is based on the intrinsic density, a concept, which requires 
additional approximations. In particular the self-energy used in DFT theory does not depend on 
the energy. Landau-Migdal theory on the other side restrains itself from calculating ground state 
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properties, but it goes far beyond the mean field approach and takes into account couphngs to 
complex configurations. There are also similarities for these two methods. Both are based on a 
single-particle description, i.e. on the motion of independent particles. In density functional theory 
an average field in introduced as a vehicle in order to take into account shell effects. However, the 
single-particle energies themselves are not observables in the strict sense. Landau-Migdal theory 
uses quasiparticles as the exact eigenstates of the A ± 1-systems. In both cases the self energies 
are given as the first derivatives of the total energy with respect to the density and the effective 
interactions between the particles are the second derivatives of this quantity. 

The combination of CDFT and TFFS described in this manuscript starts from the covariant 
density functional. No further parameters are needed. This functional is used to describe the 
ground state properties and the parameters of the functional are adjusted to experimental data 
of ground states of several nuclei. On this level the self-energy does not depend on the energy 
and Landau-Migdal theory is used to introduce an energy dependence with a partiele-vibrational 
scheme. The properties of the phonons needed for the calculation of the energy-dependent part 
of the self-energy are phonon energies and phonon-nucleon vertices. They are calculated with the 
static effective interaction obtained as the second derivative of the density functional. No additional 
parameters are needed. The essential equation of the Landau-Migdal theory is the response equa- 
tion. The effective interaction to calculate the full response is the derivative of the self-energy with 
respect to the density and this means that one obtains in addition to the static interaction resulting 
from the energy-independent part an induced interaction resulting from the energy-dependent part 
of the self-energy. Of course there would be double counting, because many of the correlations 
induced by the coupling of virtual phonons have also contributions at the ground state energy. 
Therefore a subtraction method is introduced, which removes from the induced interaction at fi- 
nite energy its value at zero energy, i.e. after this subtraction, the induced interaction vanishes at 
the ground state and takes into account only its energy dependence. Therefore this subtraction 
procedure guarantees that one does not need to readjust the parameters of the density functional, 
because the effects of partiele-vibrational coupling vanish at the ground state of the even-even 
system, where the parameters are adjusted. 

We have discussed several applications of this method, as the fragmentation of single-particle 
energies in odd-mass nuclei in the vicinity of a double magic configuration. Close to the Fermi 
surface there is always a dominant pole with a reduced single-particle strength and many other 
poles with rather small strength. The dominant pole is shifted in the direction of the Fermi surface, 
i.e. the level density at the Fermi surface is increased. The effective mass derived from this level 
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density in the Pb region is considerably increased, but there is still room for an additional energy 
dependence of the self-energy not taken into account by the coupling to surface vibrations. In 
addition we discussed several solutions of the response equations for nuclei in an external field. 
This allows to calculate the strength functions with respect to an external operator and the photo 
absorption cross sections in the correlated system. If one takes into account only the static part 
of the interaction one finds the usual RPA or QRPA results of time-dependent density functional 
theory, which reproduces the position of the resonances rather well, but it cannot account for the 
width that has its origin in the coupling to more complicated configurations. The energy-dependent 
part of the interaction includes this coupling and therefore it induces again a fragmentation of the 
rather sharp resonance peaks in RPA of QRPA over many complex configurations. Because of the 
subtraction procedure the position of the resonances is not changed very much, but the width is 
considerably increased, in excellent agreement with experimental data. 

Of course there is room for additional improvements. So far, for numerical simplicity, exist 
only applications of this theory with the parameter set NL3, which has no density dependence 
in the is isovector channel. At present there exist more modern parameter sets with the density 



dependent meson exchange 



85| which give already on the mean field level improved results for 



characteristic properties such as the neutron skin of neutron-rich nuclei or the density dependence of 
the symmetry energy. They should also be implemented in the theoretical investigations of the type 
discussed in this article. A further improvement can be achieved by using a more realistic pairing 
force. So far there are only investigations available with a monopole pairing force. It is therefore 
highly desirable to implement in the pairing channel a density dependent zero range force or the 
finite range Gogny force. All the calculations presented here have been done in a discrete basis, i.e. 
in the spectral representation of the response equation and therefore the coupling to the continuum 
is not taken into account properly so far. This might have a strong influence on application to 
light nuclei and therefore relativistic continuum RPA should be extended to relativistic continuum 
QRPA and relativistic continuum QTBA. Finally, so far ground state correlations have been taken 
into account only on the RQRPA level. In the non-relativistic case there exist investigations going 
beyond this limitations. They should be extended also to the relativistic calculations. 
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